# Presented By ZhaoEnxing
# time: 2021/10/20 10:12
# 在制药业中,药品重量的方差是很关键的。
# 对于某种特定的药物,18个单位的样本得出的样本方差S=0.36(克)。
# 构造该药物方差的90%置信区间概率。
import pandas as pd
import numpy as np


from scipy import stats
# 已知条件
confidence = 0.90
sample_n = 18
sample_var = 0.36
# 计算卡方分布（自由度n-1）90%置信区间的上下限
chi2_lower, chi2_higher = stats.chi2.interval(confidence, df = sample_n -1)
# 计算((n-1)*样本方差)的值
sample_df_var = (sample_n - 1)*sample_var
# 计算得到总体方差的置信区间估计
CI_weight_var = (sample_df_var/chi2_higher,sample_df_var/chi2_lower)

print(CI_weight_var)
